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O . Abstract 

In ferromagnetic spin models above the critical temperature (T > Tcr) the par- 
^ . tition function zeros accumulate at complex values of the magnetic field (He) with 

. a universal behavior for the density of zeros p{H) ~ |ff — He\'^ ■ The critical ex- 

■ ponent a is believed to be universal at each space dimension and it is related to 
the magnetic scaling exponent yh via a = (d — yhl/Vh- In two dimensions we have 

^ ■ yh = 12/5 (cj = —1/6) while yh = 2(a = —1/2) in d = 1. For the one dimen- 

■ sional Blume-Capel and Blume-Emery-Griffiths models we show here, for different 
temperatures, that a new value yh = 3 {a = —2/3) can emerge if we have a triple 

^ ■ degeneracy of the transfer matrix eigenvalues. 
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1 Introduction 

In [1] C. N. Yang and T.D. Lee have initiated a method to study phase transitions based on 
partition function zeros. Since then, several apphcations have appeared in equihbrium and 
more recently in non-equilibrium statistical mechanics, see [2] for a review and references. 
It is expected that with the increase of computer power new applications will appear 
specially where analytic methods can hardly be used. In many statistical models the so 
called Yang-Lee zeros (YLZ ) lie on arcs on the complex u-plane, u = exp(— where H 
is an external magnetic field and (3 = 1/KbT. The zeros at endpoints of the arcs tend to 
pinch the positive real axis ^{u) > at the critical value of the magnetic field for T < T^r 
while for T > Tcr they accumulate at complex values ue = exp(— /^if^) which are called 
Yang-Lee edge singularities (YLES ). The linear density of zeros has a universal power 
like divergence p{u) ~ \u — ue\'^ in the neighborhood of the YLES . The critical properties 
of the YLES are described by a non-unitary tip^ field theory [3] which corresponds in two 
dimensions |1] to a {p, q) = (2, 1) non-unitary conformal minimal model leading to the 
prediction a = —1/6. This is in agreement with numerical results for the ferromagnetic 
Ising model, see e.g. [5[ [6], and experimental magnetization data [?]. Even in the anti- 
ferromagnetic Ising model where ue & M- numerical results [8] also lead to a = —1/6. 
Analogously, for the ferromagnetic and anti- ferromagnetic 5 = 1/2 Ising model in ci = 1 
one derives exactly cr = —1/2. In the ferromagnetic case, the value cr = —1/2 has been 
found in other d = 1 models like the n-vector chain [9] and the q-state Potts model with 
g > 1 and < g < 1 [lO]. Furthermore, in pT], by assuming that the transfer matrix 
eigenvalues are real in the gap of zeros and that the magnetic field is pure imaginary, one 
obtains cr = —1/2 for the S = 1, 3/2 Ising model and also for the Blume-Emery-Griffiths 
(BEG) model investigated here. The same assumptions of [TT] have been made in the 
proof of cr = — 1 /2 given in [121 E] • More recently [13] another proof which holds without 
those assumptions has appeared in the Blume-Capel model. We show here that the latter 
proof contains also some hypotheses which are not satisfied in part of the parameter space 
of the BEG model where a new critical behavior appears. 



2 General setup and analytic results 

Our starting point is the couple of finite size scaling (FSS) relations, see p3l [15] . 
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p{L) = C2Ly^-^ + 



(2) 



where Ci , C2 are constants, independent of the size of the lattice L. The Yang-Lee zero 
Ui{L) is the closest zero to the edge singularity ue = Mi(oo). The density of zeros p(L) is 
calculated at u = ui{L). In practice we use 

N\ui{L} - U2{L}\ 

where U2{L) is the second closest zero to m^; and is the total number of spins (A^ = L'^). 
The dots in ([I]),© stand for corrections to scaling. Our analytic method consists of 
deriving yh from an expansion for large number of spins similar to ([1]) obtained directly 
from the transfer matrix solution. We also calculate the YLZ numerically for finite number 
of spins and check the validity of ([T]),([2]) which furnish numerical estimates for yh. 
The partition function for the BEG model is given by [I^ 



Zr, = J2expplJJ2s^S, + kJ^S^S] + ^ [H S, + A (1 - S^)] \ , (4) 
{s.} ( (ij) {ij} i=i J 

where Si = 0, ±1 and the sum is over nearest-neighbor sites. We assume ferro- 

magnetic Ising coupling J > but K , A can have any sign. All the couplings J , K , A 
are real while the magnetic field H may be complex. The temperature is defined with 
respect to the ferromagnetic coupling via the compact parameter c = exp(— /?J), where 
< T < 00 corresponds toO<c<l. In d = 1, using periodic boundary conditions, 
Si = Si+N, the model can be easily solved via transfer matrix : 

Zn = c-^ [Af + A^ + Af ] , (5) 
where the eigenvalues Ai , A2 , A3 are the solutions of the cubic equation 

A^ - asA^ + aiA - ao = , (6) 

with coefficients 



ao = hi{l-c^)[h{l + c^)-2c] 

ai = h^{l-c^) + Ai{h-c) (7) 
a2 = X + Ah 
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We have defined for convenience 



b = exp{pK) , X = exp{pA) , x = xc = exp p{A - J) (8) 

A = u + l/u = 2 cosh{f3H) (9) 

Tlie quantity u = exp{—PH) plays tlie role of fugacity in a lattice gas pLj. The partition 
function Z]\f is proportional to a polynomial of degree 2N in the fugacity. Therefore, 
all relevant information about Zn is contained in its zeros ZN{uk) = , k = 1, ■ ■ ■ , 2N. 
Due to Z2 symmetry, Z]\f{u) = Zn{1/u), half of the zeros are the inverse of the other 
half. The exact position of the YLZ can hardly be found even in one dimension. In the 
thermodynamic limit the partition function will vanish whenever two or more eigenvalues 
share the largest absolute value, see pTlfTSlfTO] . For large finite number of spins the same 
conditions are still useful, see e.g. [19], for locating the YLZ. Thus, following [13j, let us 
assume that there exists a function A^ip) = u{(p) +u^^{ip) such that when we plug it back 
in (E]) we have: 

A2 = e'^Ai (10) 
IA2I = |Ai| > IA3I (11) 

Therefore, for a large number of spins we can neglect the contribution of A3 and write 
Zn^\Xi\^ {l + e'^^). The 2N YLZ are found from 

u^iifk) = ^ [A{ifk) ± V^(<^fc)2-4] , (12) 

with (fk = {2k — , k = 1,2, ■ ■ ■ , N. In order to find ^(v^) we start from the relations 
ao = A1A2A3 , ai = A1A2 + A1A2 + A2A3 , 02 = Ai + A2 + A3. Implementing the condition 
f fTOj) and eliminating A3 we obtain: 



a, = Aia2(l + e*'^) - A?(l + e*'^ + e*2'^) (13) 
ao = Xla2e'^ - \\ [e'^ + e*^^) (14) 

Manipulating (fT3!) and (fT4|) we further eliminate Ai and find an equation for A{if\. 



Oq (1 + 2 cos v?)^ + 4 cos^ — (a^ + a^d^^ — a\a\ — 2(1 + 2 cos </)) (2 + cos </)) aociia2 = (15) 

2 
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which is the same expression obtained in [13] for the Blume-Capel model (6 = 1). At 
this point, in principle, the Yang-Lee zeros are determined by (ITB]) . once we check fill I) . 
After taking the continuum limit in the corresponding solution u{ipk) one could derive 
the density of zeros and obtain the exponent a. However, in practice, the expression 
f lT^ is a fourth degree polynomial equation for A{(p) whose solutions are cumbersome in 
general. It is not feasible to substitute them in the solutions of the cubic equation ([HD, 
which are a bit complicated too, and then check the condition (flT!) for arbitrary values of 
the parameters b,c,x. 

Fortunately, in order to compute the exponent a or yh we only need to study the 
vicinity of the YLES. The YLES is located at A{(p = 0) since at = the two, presum- 
ably, largest eigenvalues of the transfer matrix coincide (Ai = A2) which guarantees, see 
the original works [20],|21] and [11] for a recent work, the existence of phase transition. 
Therefore, only the behavior of the solutions of the equation f|T5|) about = is required, 
as already noticed in [13], remembering of course that flTT]) must hold. Nevertheless, dif- 
ferently from those authors we do not take the continuum limit. By keeping finite we 
will be able to study the large A^ expansion for the closest zero to the YLES and determine 
the exponent yh via comparison with ([1]). 

In order to scan the parameter space of the BEG model for a possible new critical 
behavior at edge singularities we have found more useful instead of using ( |T5l) to take 
consecutive derivatives of expressions (|T3l) . f|T4l) and make some mild hypotheses as shown 
next. 

First of all, from the fact that the YLES is located at A{(p = 0) it is natural to assume 
that the smallest phase among all (pk, i-e., fi = tt/N, corresponds to the closest zero to 
the YLES, namely ui{N) = u{(fi). The zero u{tt/N) is determined from A{tx/N) via 
( IT2|) . Thus, all we need is a large A^ expansion for A{tt/N). Our basic hypothesis is the 
existence of a region in the parameter space of the BEG model where we are allowed to 
expand A^i^i) in a Tayloi^] series about (^i = 0: 



A{ipi) = A(0) + ifii —1^=0 + ^3-tI^=o + • • • (16) 



dA ^jd^A^ 
^'^=° + 2\d^^ 

In order to proceed further we would like to collect information about d^A/dip"^ at <y9 = 0. 
From the first derivatives of (IT^ and fll4p we deduce: 



[5.(6 _ c) - 6 Ai] ^ = [e^^ - 1] 1^ [(2 + e^^) Ai - a^] (17) 



^Notice that due to the square root in ([T^ . when the YLES is located at A{Q) = —2 {u = —1), the 
existence of (ITBl) does not always lead to a Taylor expansion for ui{N), as assumed in [T3] 



5 



Assuming that dXi/dip at ip = is smooth enough to guarantee that the right handed 
side of (fTTI) vanishes at = we get: 



[x{b - c) - 6Ai]^^o — 1^=0 = 



dA 

Therefore if Ai 7^ x{b — c)/b we derive {dA/d{p)^^Q = 0. Since at ^9 = we have double 
degeneracy Ai = A2, back in the cubic equation ((61) it is possible to show that Ai = 
A2 = x{b — c)/b requires x = (1 — c^)b'^ /\b — c\. Thus, if x 7^ (1 — c^)6^/|6 — c| we havi^ 
{dA/ dip)^^Q = 0. Consequently, the second derivatives of (fT3|) and (IT^ furnish 



[x{b-c)~bXi]^^Q -^W=o 



(19) 



Y(Ai-A3) 



J </3=0 

Then, if there is no triple degeneracy Ai = A2 7^ A3 and x 7^ (1 — c^)6^/|6 — c| we end 
up with [dAj di^)^^^ = but [d"^ A / dip'^) 7^ 0. Back in (fT6l) and then in (|T2l) we derive 
by comparison with ([T]) the usual result yt = 2 and cr = —1/2 which generalizes for the 
BEG model the proof of o" = —1/2 given in [13] for the Blume-Capel model. On the 
other hand, if the couplings are such that we have triple degeneracy, using {dA/dip)^^^ = 
= {d"^ A / dip"^) it is easy to show that {d^ A / dip^) 7^ 0, consequently we have a new 
critical behavior with yh = 3 and a = —2/3 at the edge singularity. For completeness we 
mention the third (1 - c^)6^/|6 - c| which leads to {dA/d(p)^^Q 7^ and yt = 1, 

consequently yh — d = 0. 

Regarding the exponent yh, the point ue = —1(^4(0) = —2), which corresponds to 
I3He = ±^vr, is special as we see from (|T2l) . Due to the square root in ( fT2l) we expect a 
change from yh to yh/2 when the edge is located at ue = —1. The values yh = 1,2,3 
will be replaced by yh = 1/2, 1, 3/2 respectively. Indeed, all such special values for yh at 
Ue = —1 have been confirmed numerically, though they are not presented here with the 
exception of yh = 3/2 which appears in table 1 (third line from the bottom) for which 
a = —1/3 instead of —2/3. Even in the cases where yh = 1,1/2 our numerical results 
based on the FSS relations ([1]) and ([2]), see 

ni 

to confirm the relation a = {d — yh)/yh- It 
may be appropriate to mention at this point that there exists an analogous special point 
in the 2D Ising model for which the natural variable is u = = exp —2/3H. Namely, 
the point u = —1{/3He = ±i7r/2) is also special, see [5], there appears a change from 



^It is clear from the cubic equation ([6]) that in the special case b = c the condition Ai x(b — c)/b is 
always satisfied for any finite non vanishing temperature. 

■^For the cases where yh = 1/2, for odd number of spins, the zeros follow a more complicated pattern 
(do not lie on a continuous curve apparently) and the numerical uncertainties are larger. 
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a = —1/6 to —1/9 at u = —1. Henceforth, we concentrate on the triple degeneracy case 
(a = -2/3). 

3 Numerical results 

In this section we complement our analytic results in favor of a new critical behavior with 
a numerical study. First of all, we note that Ai = A2 = A3 = 02/8 holds if and only if : 

aia2 — 9^0 = > 3ai — 02 = (20) 

For the 5 = 1 Ising model (6 = 1 = x) there is no real temperature, not even in the 
antiferromagnetic case (c > 1), at which conditions fl20|) hold, similarly for the g = 3 
states Potts model where h = and x = u/c'. However, for the BEG model there 
are many solutions for the triple degeneracy conditions (I2UI) . In what follows we present 
numerical results for two special cases. 

3.1 BEG model with b = c{K = - J) 

For the BEG model with b = c{K = —J) the dependence of the cubic equation on the 
magnetic field simplifies which allows us to find the YLZ for a larger number of spins. In 
this case the conditions (120|) imply the fine tuning 

x = x(c) = (1 + c^) v/(l + c2)/(27(l-c2)) (21) 

and the position of the edge singularity: Ae = [2(4c^ — 5)/(l + c^)] x(c). The singularity 
is on the left-hand endpoint of an arc on the unit circle (fig. la) for < c < and 
on M_ for V2/2 < c < 1, see figures la and lb. We have checked that the YLES on the 
right-hand endpoint of the arc is of the usual type a = —1/2. 

In finding the zeros we have started from the expression for the partition function of 
the BEG model obtained by a diagrammatic expansion similar to what has been done 
in [19] for the Blume-Capel model. Following the same steps of [1^, we first define the 
generating function of partition functions of the BEG model on non-connected rings and 
then we take its logarithm which gives rise, as is well known in diagrammatic expansions, 
to the generating function of only connected diagrams (usual rings). Those steps lead us 
to = —Nc'^ — + dig"^ — aog^]}gN- The notation {f{g)}gN indicates the 

power g'^ in the Taylor expansion of f{g) about g = 0. The constant g plays the role of a 
coupling constant in a diagrammatic expansion a la Feynman (perturbation theory). The 
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above expression for Zj^ is equivalent to ([SD for arbitrary number of spins. In particular, 
the reader can easily check it for low number of spins by using the relationships between 
eigenvalues and coefficients of the cubic equation given in the text between formulae f|T2l) 
and f|T3|) . The new expression turns out to be more efficient than ([5]) from a computational 
point of view. In practice, we have used rings with Na = 210 + 10(a — 1) spins where 
1 < a < 10. The zeros have been found with the help of Mathematica software. They 
agree up to the first 30 digits with few cases where they are known exactly as, for instance, 
X = c(3 — c^)^/2 and b = c(3 — c^)/2 where, for even number of spins, the Yang-Lee zeros 
are exactly located at = -2 - 2 i{l - c^) sin (27r(l + 3k) /3N) , k = 0, ■ ■ ■ , N - 1. 

According to the discussion at the end of the previous section, we expect i/h = 3/2 
at c = V2/2 where Ae = —2, and yt = 3 elsewhere. This is approximately confirmed 
in table 1. The estimates y\l^ and y^^^ refer respectively to di-log fits of ([1]) and ([2]) with 
L = N and d = 1. We take the logarithm of the absolute value of the imaginary (real) 
part of ([T]) for < c < v^/2 (a/2/2 < c < 1) and neglect corrections to scaling in ([T]) and 
([2]). For both fits we have calculated ^-^^d the deviation from the ideal (±1) Pearson 
coefficient, but since they have the same temperature dependence we only display 
our tables. For 0.1 < c < 0.5 the results are quite homogenous but the proximity of the 
point V2/2 ^ 0.707 leads to a crossover behavior at c = 0.7 with an increase of 4 orders 
of magnitude in x^- 

Alternatively, from ([2]) with d = 1 and L = N, neglecting corrections, we obtair0: 

y^>:'\Na) = 1 + 

Some results are shown in table 2 altogether with the ^ 00 Bulirsch-Stoer (BST) 
[221 [23] extrapolation. We remark that in the BST approach there is a free parameter w 
which depends on how we approach the N ^ 00 limit, namely, yh{N) = yh{oo) + + 
+ ■ ■ ■. The BST algorithm approximates the original function yh{N) by a sequence 
of ratios of polynomials with a faster convergence than the original sequence yh{Na). 
Following [23] we can determine w in a such way that the estimated uncertainty of the 
extrapolation is minimized. By defining the uncertainty as twice the difference between 
the two approximants just before the final extrapolated result, see [23], we have tested w 
in the range 0.5 < w < 2.5 and concluded that w = 1 is the best choice. All BST results 
in this work have been calculated at w = 1. For 0.1 < c < 0.6 the extrapolations are 
remarkably close to yt = 3. The huge uncertainty at c = 0.7 (last row in the third column 
in table 2) signalizes again the crossover to yt = 3/2. Our numerical results confirm both 

''A similar formula can be deduced from ^ but it gives results even closer to yh — 3. 



In 



a+l 



In 



p{Na. 



+1) 



PiNa) 



1,2, 



(22) 



8 



finite size scaling relations ([T]) and ([2D. See fig. 2 for a plot regarding the relation ([2]). 

Before we go to the next case we can, without much numerical effort, check that 
the triple degeneracy conditions lead in fact to new results. When we choose b = c the 
equation f|T5l) becomes a cubic equation for A{ip). By further fixing x = x{c) according 
to fl?I]) and c = 0.4, for sake of comparison with the zeros in fig.l, we find three solutions 
Ai{(p) , z = 1, 2, 3 which are still complicated, but expanding about if = we derive: 



, , , 829 /29 9^203 ^ ^203 4 . 

900 V 7 400 ^ 1440 ^ ^ 

^ , , 218 /29 1 [203 3 V203 4 . 5. 

225 V 7 75 V 3 ^ 450 ^ ^ 

A3(^) = A2(-^) (25) 



(23) 
(24) 



Now we see that the symmetry ip ^ —ip oi the equation ( |T5l) mentioned in [13] does not 
necessarily lead to even solutions. The operation (f — > —ip exchange the second with the 
third solution. The first solution Ai{(p) is the usual Yang-Lee edge singularity and gives 
cr = — 1/2 while the second one leads to a = —2/3. The second solution corresponds to the 
new edge singularity appearing on the left side of the arcs in fig. la. The third solution 
is not a true singularity. We can, even before the expansion about if = 0, substitute 
Ai{ip) in the eigenvalue solutions of the cubic equation and check numerically that 
the magnitude of two eigenvalues degenerates for all Ai{ip). However, for ^3(^9) such 
magnitude is smaller than the real (non-degenerate in magnitude) eigenvalue and so it 
does not satisfy ( fTTl) . 

3.2 Blume-Capel model (b = 1) 

For the Blume-Capel model K = 0{b = 1) the triple degeneracy conditions ( l20l) lead to 
complicated formulas for both the fine tuning function x = x{c) and the position of the 
new edge Ae ■ 



x{c) = {1 - 8c + 10c' + lOc^ - 9c^ + 3(1 - c') [{P + Qf'^ +R{P + Q)"^/' ] }'^' (26) 

^ x{c) [(1 + c)(5 - 15c + 5c' + 3c^) + x^jc)] 

(l + c)(l + c2) + x2(c)(2-3c) ^ ^ 

Where P = y/c^{l + c - c'^){l - 3c - , Q = c2(5 - 36c + 95c2 - 90c=^ + 27c^) and 
R = c(9c^ — 20c' + 10c — 1). In the whole range of temperatures < c < 1 the function 
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x{c) given in fl26l) runs from +1 to +2, in fact we have approximatelju x(c) ~ 1 + c. The 
position of the zeros Ae ranges from -2 to +2 which guarantees, see that the new edge 
is always on the unit circle on the complex u-plane. However, the finite size numerical 
calculations reveal that the zeros on the unit circle correspond to a small fraction of whole 
set of zeros, see fig. 3 . Those zeros tend to form a short arc whose right-hand endpoint 
is a usual YLES with cr = — 1/2 as we have checked. The unusual behavior appears on 
the left-hand endpoint, whose results are displayed in tables 3 and 4. There is no edge 
singularity for the other zeros outside the unit circle. For technical numerical reasons we 
have used Na = 96 -|- (a — 1)6 spins with 1 < a < 10 and 0.3 < c < 0.9. Although the finite 
size results move away from i/h = 3 as T (table 3), the EST extrapolations (table 4) 
all tend to yh = 3 with increasing uncertainty as T — > 0. In table 4 we only display the 
limit cases c = 0.3 and c = 0.9. Since we have noticed that the fraction of zeros on 5*1 
decreases as T — 0, our interpretation is that we are closer to a continuum distribution of 
zeros as the temperature increases and the temperature dependence of yll'\N) is a pure 
finite size artefact. In conclusion, our numerical results are in good agreement with our 
analytic predictions. 

Analogously to what we have done at the end of last subsection, if we now assume 
b = 1, X = x{c) as given in (l26l) and fix some specific value for the temperature, we 
show that the expansions of the solutions of the equation (|T5l) around (p = give both 
a = -1/2 and a = -2/3. 



4 Conclusion 

We have derived for the Blume-Capel and BEG models, analytically and numerically, 
a new value for the magnetic scaling exponent yt = S and edge exponent a = —2/3. 
This is only possible due to a triple degeneracy of the transfer matrix eigenvalues. We 
have checked that our di-log fits, see e.g. fig.2, are consistent with the finite size scaling 
relations ([1]) and ([2]) which appear around the usual Yang-Lee edge singularity. We stress 
that the new critical behavior found here is not in conflict with previous numerical and 
analytic results in the literature. In particular, the proof of a = —1/2 given in [T^l El ITT] 
assumes that in the gap of zeros close to the YLES all eigenvalues are real. In the 
examples that we have analyzed so far we have verified that this is indeed the case for 

^Although i(c) w 1 + c it is important that x{c) ^ 1 + c since 5; = (1 — c^)fe^/|6 — c| = 1 + c for & = 1 
and, as we have discussed after formula (jlSp . at that point we may non longer have dA/dip — at ip = 0. 
In fact, we have confirmed numerically that « 1 at this point. It corresponds to a bifurcation point 
where one curve becomes two disjoint arcs of zeros, see fig. (2c) of [19j 
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the usual YLES (a = —1/2), however only the largest eigenvalue is real in the gap close 
to the new singularity (a = —2/3). Besides, although all eigenvalues are real at vanishing 
magnetic field, as one can check from ([6]) at A = 2, apparently nothing prevents the two 
eigenvalues with smallest magnitude of acquiring an imaginary part before reaching the 
edge singularity as we increase the magnetic field through imaginary values. Consequently, 
so far we have not been able to prove that the edge singularity closer to M+ is always of 
the usual type (cr = —1/2), or alternatively, that the triple degeneracy condition could be 
satisfied also by an edge singularity close to the positive real axis. Concerning the recent 
proof of cr = —1/2 given in [13], it supposes that {dPu/ d({P')^^Q ^ which is not the case 
in the examples presented here where Ai = A2 = A3 7^ x{h — c)/h. 

Finally, it is important to remark that instead of imposing the triple degeneracy condi- 
tions (|20|) and finding rather peculiar fine tuning functions fl2Tl) and fl26l) we could turn the 
argument around and claim that for some given values of the couplings h and x we can, in 
many cases, find a specific real temperature c such that conditions fl2I]]) are satisfied and 
the usual result a = —1/2 is replaced by a = —2/3. Since the number of transfer matrix 
eigenvalues is infinite for > 1, even for the S* = 1/2 Ising model, we might speculate 
that the degeneracy of more than two largest eigenvalues could play a role also in higher 
dimensional spin models and lead to a change in the exponent a at specific temperatures 
in (i > 1. 
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Fig. la Fig. lb 




Figure 1: Yang-Lee zeros for the BEG model with 6 = c on the complex u-plane {u = 
exp~^^) at the triple degeneracy condition, see (12T|) . and temperatures c = 0.4 (Fig. la), 
c = 0.9 (Fig. lb). In both figures we have 100 zeros. The closest zero to the new edge 
corresponds to uf ^ -0.98602 ± 0. 16660 i in Fig. la and ^ -0.74063 , ^ -1.35021 
in Fig. lb 



c 






2 

Vh 


lO^v^ 

J-U X(2) 


0.1-0.5 


2.997 


3.5 


2.989 


4.0 


0.6 


2.997 


3.5 


2.987 


8.3 


0.7 


2.880 


3.2 xlO^ 


2.094 


2.4x10^ 


V2/2 


1.498 


0.9 


1.494 


1.3 


0.8 


3.000 


3.3 


2.989 


2.3 


0.9 


2.997 


3.5 


2.989 


3.8 



Table 1: Data from di-log fits of formulae ([T]) and ([2]) with L = N and d = 1 for the 
A = - J (6 = c) BEG model with 210 < < 300 spins. 
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ln[p] 




In [N] 



Figure 2: Di-log fit for tlie density of zeros around tlie new edge singularity for the 
A = — J {b = c) BEG model. We have used 210 < N < 300 spins at temperature c = 0.6. 
The sohd line is Inp = 1.98724 (In A^) - 6.63547. 



Na 


1/^(0 = 0.6) 


yi'\c = 0.7) 


210 


2.98450354337 


1.96268002086 


220 


2.98540244520 


1.99853987788 


230 


2.98619971860 


2.03403890400 


240 


2.98691163425 


2.06907936446 


250 


2.98755121168 


2.10357440914 


260 


2.98812897995 


2.13744702199 


270 


2.98865353760 


2.17062929861 


280 


2.98913197003 


2.20306194085 


290 


2.98957016460 


2.23469387995 


oo 


3.000009(6) 


3(8) 



Table 2: Finite size results, obtained from (l22l) . for the magnetic scaling exponent yh for 
the {b = c) BEG model with 210 < N < 300 spins at temperatures c = 0.6 and c = 0.7. 
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Figure 3: Half (Im(M) > 0) of the Yang-Lee zeros for the Blume-Capel model (6 = 1) on the 
complex u-plane {u = exp~^^) at the triple degeneracy condition, see (126|) . temperature 
c = 0.5 and N = 150 spins. The solid line in both figures is part of the unit circle. The 
closest zero to the new edge corresponds to ~ 0.2332 + 0.9734 i as shown in detail in 
Fig.3b. 



c 






(2) 

Vh 




0.3 


3.200 


5.03 


4.612 


3.10 


0.4 


3.148 


2.52 


3.795 


2.10 


0.5 


3.119 


1.51 


3.561 


0.72 


0.6 


3.100 


1.00 


3.442 


0.36 


0.7 


3.088 


0.76 


3.370 


0.22 


0.8 


3.078 


0.58 


3.320 


0.15 


0.9 


3.071 


0.48 


3.285 


0.11 



Table 3: Data from di-log fits of formulae ([T]) and ([2]), with L = N and d = 1, for the 
Blume-Capel model with 96 < < 150 spins and temperatures 0.3 < c < 0.9. We have 
used \^e{ut{N) - mi(cx)))| in the di-log fit of (p. 
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yl'\c = o.3) 


yi'\c = 0.9) 


96 


6.65138835174 


3.36602639165 


102 


5.40615603620 


3.33723696161 


108 


4.85302919346 


3.31271060187 


114 


4.52543971437 


3.29155556642 


120 


4.30452579272 


3.27311462040 


126 


4.14377246532 


3.25689223008 


132 


4.02074088023 


3.24250742353 


138 


3.92311802198 


3.22966230906 


144 


3.84352213555 


3.21812047991 


oo 


3.001(1) 


3.00000(8) 



Table 4: Finite size results, obtained from (l22l) . for the magnetic scaling exponent yh for 
the Blume-Capel model {b = 1) with 96 < < 150 spins at temperatures c = 0.3 and 
c = 0.9. 
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